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Abstract 

A solution tp to Schrodinger's equation needs some degree of regularity 
in order to allow the construction of a Bohmian mechanics from the inte- 
gral curves of the velocity field HQ (\/ip/mil>) . In the case of one specific 
non-differentiable weak solution ^ we show how Bohmian trajectories can 
be obtained for ^ from the trajectories of a sequence \]/ n — > ty. (For any 
real t the sequence (t, •) converges strongly.) The limiting trajectories 
no longer need to be differentiable. This suggests a way how Bohmian 
mechanics might work for arbitrary initial vectors \I/ in the Hilbert space 
on which the Schrodinger evolution \1/ i-> e~ tht} i> acts. 

1 Introduction 

Quantum mechanis often is praised as a theory which unifies classical mechanics 
and classical wave theory. Quanta are said to behave either as particles or waves, 
depending on the type of experiment they are subjected. But where in the 
standard formalism can the particles of the interpretive talk be found? Perhaps 
only to some degree in the reduction postulate applied to position measurements. 
In reaction to this unsatisfactory state of affairs, Bohmian mechanics introduces 
a mathematically precise particle concept into quantum mechanical theory. The 
fuzzy wave functions are supplemented by sharp particle world lines. Through 
this additional structure some quantum phenoma like the double slit experiment 
have lost their mystery. 

Clearly the additional structure of particle world lines brings along its own 
mathematical problems. Ordinary differential equations are generated from 
solutions of partial differential equations. A mathematically convincing general 
treatment so far has been given for a certain type of wave functions which do 
not exhaust all possible quantum mechanical situations. Exactly this fact has 
led some workers to doubt that a Bohmian mechanics exists for all initial states 

of a Schrodinger evolution t h-» e~ lht ^Q. We shall show on one specific case 
of a counter example ^f now the problem might be resolved in general. We 
approximate the state for which the Bohmian velocity field does not exist, 



by states which do have one. Their integral curves turn out to converge to limit 
curves which can be taken to constitute the Bohmian mechanics of the state 
unamenable to Bohmian mechanics on first sight. 



2 Bohmian evolution for ip G C 2 

Let ip : R x R s — > C be twice continuously differentiable, i.e. ip S C 2 (R x R s ) , 
and let ip obey Schrodinger's partial differential equation 

h 2 

iHd t ip(t,x) = -—Aip(t,x) + V (x) tp (t, x) (1) 

with V : R s -> R being smooth, i.e., V € C°° {R s ) . From ip, which is called 
a classical solution of Schrodinger's equation, a deterministic time evolution 
x ^ lx (*) °f certain points a; € R s can be derived: If there exists a unique 
maximal solution r y x : I x — > R s to the implicit first order system of ordinary 
differential equations 

M*,7(i))7(*)=^(*.7(t)) ( 2 ) 

with the initial condition 7 (0) — x, one takes 7^ as the evolution of x. Here 
^ :lx R s -> R and j^RxiM R s with 



Pii (*> a; ) = IV'C*)^)! 2 and j,/, (i, m) = — 9 tp (t,x)V x rp (t,x) (3) 

m 



obey the continuity equation dtp^, (t,x) = — diyy',/, (i, x) for all (t, x) e R x I s . 
From now on we shall drop the index if) from and j^. 

For certain solutions^] tp the curves j x can be shown to exist on a maximal 
domain I x = R for all x G M s : If ip has no zeros, then the velocity field 
v = j/p is a C 1 -vector field, w then obeys a local Lipschitz condition such 
that the maximal solutions are unique. If in addition there exist continuous 
nonnegative real functions a, /3 with \v (t, x)\ < a (t) \x\ + /3 (t) then all maximal 
solutions equation ^ are defined on R and the general solution 

$ : U xeR . ^ x M RS with $ (*' x ) = (*) 

extends to all of R x R s . (Thm 2.5.6, ref. [T]) Due to the uniqueness of maximal 
solutions the map $ (t, •) : R s — > R s is a bijection for all t € R. It obeya^] 

p(t,x)d s x = \ p(0,x)d s x (4) 
#(t,n) in 



1 The simplest explicitly solvable example is provided by the plane wave solution rjj (t, x) = 
e — i\k\ t+ik-x ^ Bohmian evolution <E> obeys <& (t, x) = tk. Another explicitly solveable case 
is given by a Gaussian free wave packet. 

2 Here <S> (t, f2) = {* (t, x)\xeQ}. 
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for all t € K and for all open subsets fl C M s with sufficiently smooth boundary 
such that the integral theorem of Gauss can be applied to the space time vector 
field (p,j) on the domain Ut'g(o t) ^ (*'> ^) • 

These undisputed mathematical facts have instigated Bohm's amendment 
of equation ([!]) in order to explain the fact that macroscopic bodies usually are 
localized much stricter than their wave functions suggest. 

In Bohm's completion of nonrelativistic quantum mechanics it is assumed 
that any closed system has at any time, in addition to its wave function, a 
position in its configuration space and that this position evolves according to 
the general solution $ induced by the wave function. One says that the position 
is guided by tp since $ is completely determined by tp (and no other forces than 
the ones induced by -0 ar e allowed to act on the position). More specifically, 
7 X is assumed to give the position evolution for an isolated particle with wave 
function ij) (0, •) and position x - both at time t = 0. 

As is common in standard quatum mechanics, -0 (0, •) is supposed to obey 

\ip(0,x)\ 2 d s x = 1. 

The nonnegative density p (0, •) is interpreted as the probability density of the 
position which the particle has at time t = 0. Since an initial position x is 
assumed to evolve into "f x (t) , the position probability density at time t is then, 
due to equation Q, given by p(t,-). In particular, Bohm's completion gives 
the position probabilities among all the other spectral probability measures 
a fundamental status, since the empirical meaning of the other ones, as for 
instance momentum probabilities, all are deduced from position probabilities. 

There are classical solutions of Schrodinger's equation, whose general so- 
lution $ does not extend to all of K x R s . An obstruction to do so can be 
posed by the zeros of ip. In the neigbourhood of such zero the velocity field 
v = j/p may be unbounded and v then lacks a continuous extension into the 
zero. As an example consider a time wave function ip (0, •) : R 2 — > C, for 
which ip (0, x, y) = x 2 + iy within a neighbourhood U of its zero {x, y) — (0, 0) . 
Within U for the velocity field follows 

m ! , . d x ip{0,x,y) 2xy 2 
-v {0,x,y) = 3?- 



h y ' ' y/ ip(0,x,y) x 4 + y 4 ' 

Hence for < \<fi\ < ir/2 we have v 1 (r cos <fi, r sin (j>) — > — oo for r — > with <\> 
fixed. Thus the implicit Bohmian evolution equation ([2]) is singular in a zero 
of the wave function whenever the velocity field does not have a continuous 
extension into it. As a consequence the evolution j x of such a zero x is not 
defined by equation pj). 

As a related phenomenon there are solutions to equation |2]) which begin or 
end at a finite time because they terminate at a zero of if). A nice example {3| 
for this to happen provide the zeros of the harmonic oscillator wave function 
ip : M 2 -> C with 

i){t,x) = e-^ (l + e~ 2 "(l-2x 2 )). 
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E.g., the points \x\ = 1 are zeros of ip(t,-) at the times t 6 ttZ. They are 
singularities of v since 

t->o ip(t,±l) 

Note however that S ^ffi = for a: ^ ±1. 

There are more challenges to Bohmian mechanics. The notion of distri- 
butional solutions to a partial differential equation like raises the question 
whether these solutions support a kind of Bohmian particle motion like the 
classical solutions do. After all quantum mechanics employs such distributional 
solutions. 



3 Bohmian evolution for G C£° 

In standard quantum mechanics the classical solutions, i.e. the C 2 -solutions of 
equation ([!]), do not represent all physically possible situations. Rather a more 
general quantum mechanical evolution is abstracted from equation 0. It is 
given by the socalled weak solutions 

*o H' * t = e- iht ^ for all * G L 2 (M s ) 

with hh being a self-adjoint, ususally unbounded hamiltonian corresponding to 
equation The domain Dh of h does not comprise all of L 2 (R s ) , yet it 

is dense in L 2 (M. s ) . Since h is self-adjoint, the exponential e~ lht has a unique 
continuous extension to L 2 (M. s ) . This unitary evolution operator e~ lht stabilizes 
the domain of h as a dense subspace of L 2 (M. s ) . Thus if and only if an initial 
vector * belongs to Dh, equation generalizes to 



lim 

6-S-O 



I 



(5) 



for all i£R. For ^ Dh equation ^ does not hold for any time. 

Yet the construction of Bohmian trajectories needs much more than the 
evolution ^ &t within L 2 (R s ) , since the elements of L 2 (M. s ) are equivalence 
classes [/] of functions / e C 2 (R s ) . It rather needs a trajectory of functions 
instead of a trajectory of equivalence classes of square-integrable functions. If 
there exists a function ifeC'flx R s ) such that * (t, •) = [ip (i, •)] holds for all 
tel, then ip is unique and the Bohmian equation of motion ^ can be derived 
from the evolution vJ/ i-> \f> t through ip. When does there exist such if>1 

Due to Kato's theorem (e.g. Thm X.15 of ref. [1]) the Schrodinger hamil- 
tonians h, corresponding to potentials V from a much wider class than just 
C°° (M. s : M) , have the same domain as the free hamiltonian —A, namely the 
Sobolev space W 2 (R s ) . This is the space of all those $€i 2 (R s ) which have 
all of their distributional partial derivatives up to second order being regular 
distributions belonging to L 2 (R s ) . Since Dh is stabilized by the evolution e~ tht , 
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for any ^ £ Dh there exists for any tela function ip(t,-) G W 2 (M s ) such 
that 

e- iht * = [><p (V)]. (6) 

However, for this family of time parametrized functions ij)(t,-) the Bohmian 
equation of motion in general does not make sense since ip (t, •) need not be 
differentiable in the classical sense0 

Therefore some stronger restriction of initial data than ^>o 6 Dh is needed 
in order to supply the state evolution ^ e~ lhtl $o with Bohm's amendment. 
For a restricted set of initial states (x, an d for a fairly large class of static 
potentials a Bohmian evolution has indeed been constructed in [5] and [B] . There 
it is shown that for any "Jo g C\ n eN^ hn =: there exists 

• a (time independent) subset f2 C R s 

• for any t a square-integrable function ip (t, •) 

such that the restriction of ip (t, •) to Q belongs to C°° (Q) and equation (|6| 
holds. The set Q is obtained by removing from K s first the points where the 
potential function V is not C°°, second the zeros of tp (0, •) , and third those 
points x for which the maximal solution j x does not have all of M. as its domain. 
Surprizingly, f2 is still sufficiently large, since 

f \ip{0,x)\ 2 d s x= 1. 
Jn 

On this reduced set f2 of initial conditions a Bohmian evolution $ : M x il —> W 
can be constructed. Thus if V&o G C£° and if the initial position x is distributed 
within K s with probability density (0, -)| 2 then the global Bohmian evolution 
j x of x exists with probability 1. 

4 Bohmian evolution for f ( G L 2 \ 

How about initial conditions V^o G £ 2 (IR S )\C^ ? Can the equation of motion ^ 
still be associated with Hall has devised a specific counterexample 4 Dh 
which leads to a wave function tp which at certain times is nowhere differentiable 
with respect to x and thus renders impossible the formation of the velocity 
field v. Therefore it has been brought forward that the Bohmian amendment of 
standard quantum mechanics is "formally incomplete" and it has been claimed 
that the problem is unlikely to be resolved. [7] 

A promising way to tackle the problem is to succesively approximate the 
initial condition Vf^ 4- Dh by a strongly convergent sequence of vectors (^q) 

3 Only fors=l Sobolev's lemma (Thm IX.24 in Vol 2 of ref. [J) says that [ip (t, •)] has a C 1 
representative within C 2 (M) . From such a C 1 representative ip (t, ■) the current j follows as a 
continuous vector field and a continuous velocity field v can be derived outside the zeros of tp. 
However, v does not need to obey the local Lipschitz condition implying the local uniquenss 
of its integral curves. 
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C^°. For each of the vectors a Bohmian evolution <£>„ exists. We do not 
know whether it has actually been either disproven or proven that the sequence 
of evolutions does converge to a limit $ and that the limit depends on the chosen 
sequence \&q — > <f . 

Here we shall explore this question within the simplified setting of a spa- 
tially one dimensional example. We will make use of an equation for r y x which 
has already been pointed out in [5] and which does not rely on the differentia- 
bility of j. In this case equation Q can be generalized in order to determine 
a nondifferentiable Bohmian trajectory j x by choosing Q = (—00, x) in Q as 
follows. 

Consider first the case of a C 2 -solution of equation generating a gen- 
eral solution $ : R x R -> t of the Bohmian equation of motion Q. Since 
because of their uniqueness the maximal solutions do not intersect, we have 
$ (t, (—00, x)) — (—00, $ (t, x)) = (—00, j x (£)) . From this it follows by means 
of equation Q that 

p{t,y)dy = f p(0,y)dy. (7) 



As a check we may take the derivative of equation ^ with respect to t. This 
yields 

P(t,l x (t))i x (t)+ d t p(t,y)dy = 0. 



Making use of local probability conservation d t p = —d x j we recover by partial 
integration equation 

Now observe that the equation ^ for r y x (t) is meaningful not only when tp is 
a square integrable C 2 -solution of equation (fiT) but also if i)j (t, •) is an arbitrary 
representative of with arbitrary € L 2 (K) . In order to make this explicit 
let E x : L 2 (R) — > L 2 (R) with x e R denote the spectral family of the position 
operator. For the orthogonal projection E x holds 



i.E x <p) (y) 



tp (y) for y < x 
otherwise 



The expectation value ('I', E x ty) of E x with a unit vector 4/£i 2 (R) thus 
yields the cumulative distribution function of the position probability given by 

If we define F : R 2 ->• [0, 1] through F (t, x) = (* t , , then equation Q 

is equivalent to 

F(t, Jx (t))=F(0,x). (8) 

Thus, the graph {(t, j x (t)) \t g R} of a trajectory is a subset of the level set of F 
which contains the point (0, x) . If is a sequence in L 2 (R s ) which converges 
to * then 



lim F n (t,x) = lim l^ ni e lht E x e- rht ^ n ) = lim ll^e-^^, 
= ||£ x e- iht f|| 2 = F(t^) 
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because e~ lht ,E x , and ||-|| are continuous mappings. 

Note that for any t £ K the function F(t,-) : K — » [0,1] is continu- 
ous and monotonically increasing. Furthermore linx^-oo F (t, x) = and 
linx^oo F (t,x) = 1. The monotonicity is a strict one if tp (t,-) does not vanish 
on any interval. Thus for any (£, x) S R 2 there exists at least one j x (t) £ M 
such that equation ^ holds. (For those values t for which Fit,-) is strictly 
increasing, there exists exactly one j x (t) £ K such that equation Q holds.) 
The function F cannot be constant in an open neighbourhood of some point 
(t, x) if the hamiltonian is bounded from below. Thus for any ieK, for which 
there does not exist a neighbourhood on which F (0, •) is constant, we now define 
7^ : K — > K to be the unique continuous mapping for which 

F(t, lx (t))=F(0,x). 

Note that j x : K — » K is continuous, yet need not be diffcrentiable. 



5 Hall's counter example 

Let us now illustrate this construction of not necessarily differentiable Bohmian 
trajectories by means of a solution t H> ty t 4- Dh of the Schrodinger equation 
([5| describing a particle confined to a finite interval on which the potential 
V vanishes. This solution has been used by Hall [7j as a counter example 
to Bohmian mechanics. Similar ones have been used in order to illustrate an 
"irregular" decay law t H > \ (^/q, ty t ) | 2 . 8 Both works have made extensive use 
of Berry's earlier results concerning this type of wave functions. [S] 

The (reduced) classical Schrodinger equation corresponding to the quantum 
dynamics is 

id t ip(t,x) = -d%il>{t,x) (9) 

for all (t, x) £ R x [0, w] together with the homogeneous Dirichlet boundary 
condition if; (t, 0) = -0 (t, ir) = for all t £ M. The corresponding hamiltonian's 
domain Dh is the set of all those "J £ 1? (0, 7r) which have an absolutely contin- 
uous representative ip vanishing at and tt and whose distributional derivatives 
up to second order belong to L 2 (0, ir) . As an initial condition we choose the 
equivalence class of the function 

tjj (0, x) = 1 / Vtt for all x £ [0, n] . 

Since within the class vp = [ip (0, •)] there does not exist an absolutely contin- 
uous function vanishing at and tt the equivalence class "Jo does not belong to 
Dh- As a consequence for any t the vector ^> t = e~" ! '*\E r Q does not belong to Dh- 
This in turn implies that ty t does not have a representative within the class of 
C 2 -functions on [0, tt] with vanishing boundary values. 

The hamiltonian h is self-adjoint. An orthonormal basis formed by eigen- 
vectors of h is represented by the functions Uk with 

[2 

Uk {%) = \ — sin (kx) for < x < 7r and k £ N. 
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For n e N the C°°-function ip n : K 2 -)• C with 

4 ™ e - l (2fe+i) 2 t 
V„ (t, x) = — = ^ — — sin [(2fc + 1) x] 

is a classical solution to the Schrodinger equation ^ on R 2 and fulfills homo- 
geneous Dirichlet boundary conditions at x — and x — it. Furthermore ip n 
is periodic not only in x but also in t with period 2tt. More precisely ip (t, •) is 
an odd trigonometric polynomial of degree 2n + 1 for any ( £ R. In addition 
ip n (t, •) also is even with respect to reflection at tt/2, i.e., it holds 



for all x e R. The functions (-,x) are trigonometric polynomials of degree 
(2n + l) 2 . 

As is well known, the sequence (ip n (0, -)) n eN converges pointwise on R. Its 
limit is the odd, piecewise constant 27r-periodic function cr (0, •) with 



lim Vnipn (0; x) = a/ 7 '" ' (0, a;) 



1 for < x < 7r 
for x € {0, tt} 



it (0, •) is discontinuous at x e tt ■ Z. For any t £ R the sequence (ip n (t, -)) neN 
converges pointwise on R to a function ip (£, •) . For rational £/7r this function 
is piecewise constant. [9] However for irrational t/n the real- and imaginary 
parts of ip (t, •) restricted to any open real interval have a graph with nonintcgcr 
dimension. [5] Thus for irrational t/w the function ip (t, •) is non-differentiablc 
on any real interval. As an illustration we give in Figure [T] the graph of 



x H> ^i^/inp 



o()0 



7T 

, TTX 
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for < x < 1/2 together with the partial sum over k <G {501, . . . 750} visible as 
the small noisy signal along the abscissa 



x i — y 



SR^t? (-0750 (V V^, ttx) ~ ^500 yH 12, ttx) ) 



Similarly, for given X € (0, 7r) the mapping t ip (t, x) does not belong to 
the set of piecewise C 1 -functions on [0, 2tt] . This can be seen as follows. First 
note that for given x the 27r-periodic function ip (•, x) has the Fourier expansion 

oo 

ip (t, x) = c n e~ mt where 
fe=i 

/ ^ 2FPT sin K 2/s + 1) ^ for " = (2fe + l) 2 with fc € N 
I otherwise 
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-0.4 1 1 1 1 1 1 1 1 1 1 

0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



X 

Figure 1: Real part of ip 500 at a fixed time 



Assume now that ip(-, x) is piecewise C . Then, according to a well known 
property of Fourier coefficients, there exists a positive real constant C such that 
n\c n \ < C for all n £ N. This implies that 



(2k + 1) |sin [(2k + l)x]\<C for all k £ N 



(10) 



with the positive constant C — tt^/ttC/A. However, for x 7T-Z there exists some 
real constant e > such that the set {k £ N : |sin [(2k + 1) x]\ > e} contains 
infinitely many elements. Thus for x ^ it ■ Z the estimate ( 10 ) cannot hold and 
therefore the function t\-¥ ip(t, x) cannot be piecewise C 1 on [0, tt] . 
In Figure [2] we plot the time dependence 



t I — y 3*^^15 (Trt.Tr/ 2 ) = 



A 15 

TT ^ 2k 

k=0 



-1)* 



COS 



{(2k + ifwt 



for < t < 1/2 together with the partial sum 



(noisy signal along abscissa). 

The restriction of the limit a (0, •) to [0, 7r 



V> (0, •) does. Thus lim n _ 



i>n (0, 



represents the same L 2 element as 
= 0, when ?/>„ denotes the restric- 



tion of ■0 n to K x [0, 7r] . Correspondingly p n (0, •) = \ip n (0, •) | 2 converges towards 
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Figure 2: Real part of ip 15 at x = ir/2 



the density of the equipartition on [0, tt] . Additionally, due to the continuity of 



the evolution operator e %ht , the sequence of equivalence classes 
approximates the L 2 vector ^> t = e~ lhtl $o, i.e., for all t holds 



i>n it,') 



eCg> 



lim 



= 0. 



Since also E x : I? (K) — > L 2 (R) is continuous, the time dependent cumu- 
lative position distribution function F : M x [0, 7r] — > [0,1] with F(t,x) := 
(y t ,E x y t ) obeys 



F(t,:c) = lim ( V n (*,-) 



lim / \tp n (t,y)\ dy. 

n— f oo 



The level lines of the functions F n : K x [0,tt] — > [0,1] with F n (t,x) := 
Jo IV'n (*j 2/)| 2 thus converge to the continuous level lines of F. 

Figure [3] shows some level lines of F n for n = 5, 10, 20 starting off at equal 
positions at f = 0. The level lines inherit the period 7r/4 of F (•, x) , which has 
this periode since the frequencies appearing in the even function \ip n (-,x)\ are 
0,8,16,... 

Figure [4] shows the case n = 1000. Increasing n from 20 to 1000 hardly alters 
the level lines. 
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n = 5 n = 10 n = 20 space 

Figure 3: Level lines of F n for n = 5, 10, 20 




n = 1000 space 

Figure 4: Level lines of Fiooo 
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